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(^f) , This paper presents some relations for orthonormal bases in the Minkowski space 

and isotropic tetrads constructed from the vectors of these bases. 
^3 . As an example of an application of the obtained formulae, in particular recursion 



relations, a new method is proposed to calculate traces of Dirac 7-matrices in the 
Minkowski space. Compared to the classical 
more compact expressions for the traces. Spc 
as a simple yet efficient computer algorithm. 



Minkowski space. Compared to the classical algorithms, the new method results in 
more compact expressions for the traces. Specifically, it may be easily implemented 



1 Introduction 



The problems arising when the traces of a large number of Dirac 7-matrices should be 
calculated are well known. For ten matrices and more the expressions for the traces 
become extremely cumbersome. 

For example, the classical algorithm based on the application of the algebra of 7- 
matrices 

+ (1) 

for the trace of ten 7-matrices results in the expression consisting of 945 terms (e.g. such 
algorithm is used in REDUCE [T]). 

In four-dimensional space-time somewhat shorter results provides the method using 
the Chisholm identities [2] and Kahane algorithm |3] that is based on the reduction formula 

y*YY = g ^Y - + g up Y + ie^ px ~f 5 ~f X . (2) 

This method, implemented first in SCHOONSCHIP |I] and later in FORM |5., gives for 
the trace of ten 7-matrices the following expression consisting of 693 terms: 

0x02^3^4 &5&6^7^8^'9^1o) — Tr(ttifl2^3^4^5'^6^'7^8'^9^'lo) 
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The proposed new method results in more compact expressions (e.g. consisting of 
32 terms for the trace of ten 7-matrices) than the classical methods. In particular, this 
method may be easily implemented as a simple yet efficient computer algorithm. 

The method is based on the properties of orthonormal bases in the Minkowski space 
and isotropic tetrads (see e.g. [Hj) constructed from the vectors of these bases. 

In this paper we use the Feynman metrics: 



/i = 0,1,2,3, a M = (ao,a), = (ao, —a), (ab) 
and a sign of the Levi-Civita tensor is determined as £0123 



+1; 



ao&o — ob , 
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Together with the regular Gram determinants, we use the generalized ones. For ex- 
ample: 
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2 Some relations for isotropic tetrads in the Minkowski 
space 

Let vectors Iq, h, 1%, h form an orthonormal basis in the Minkowski space, i.e. 

{IaIb)=9ab, {A, B = 0,1, 2, 3) 

(index of the basis vector is denoted by the capital letters). 

In four-dimensional space both the regular and generalized fifth-order Gram determi- 
nants equal zero, therefore 
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Whence we obtain the expression for the metric tensor: 

9 — L o L o ~ nn ~ hh ~ l 3 L 3 ■ 
Furthermore, for the orthonormal basis we have 

q ( k k k k 
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(6) 
(7) 



The "+" and "— " signs correspond to the "right-hand" and "left-hand" quadruples of the 
vectors. 

For defmiteness in the following we consider the "right-hand" basis. 



kkkk 



Tr (lohhh) + 7 5 Tr (y%khh) + Tr (a^kkU 
Tr (jHohkh) = z 7 5 W^i^a = H 5 , 



which gives us additional relations for the basis vectors, for example: 

k = irfkhh , k = il%hh , k = -rfhkk , h = h%hk 
Using (|9J). we derive the expressions which we need in the following: 



(8) 



(9) 



kh = -ij 5 kh , 


(10) 


kk = rfkk , 


(11) 


kh = -ij 5 hh , 


(12) 


kh = h%k , 


(13) 


kh = -ij 5 hh , 


(14) 


kh = il%k ■ 


(15) 



Let us introduce the isotropic tetrad 

k + k 



k — k 



(q+q- 



(16) 



e+ = ^ , e_ = , (e+e_) = -1 , (e T ) = e± 

(9+e+) = (g+e_) = (g-e+) = (g_e_) = . 
Then (0} takes the form 

III! U V i U V U V U V 

g v = q ^ q _ + q ^ q+ _ e ^ e _ _ e ^ e+ j 

from which follow the expressions below: 

(ab) = (aq + )(bq-) + (aq^)(bq + ) - (ae + )(6e_) - (ae_)(6e + ) 

7 M = Q+Q- + Q-Q+ - e +e- - e^e + , 
a = (aq + )q^ + (ag_)g + — (ae + )e_ — (ae_)e+ . 
From (fTOj) - (p"5j) we have 

q+e- = ^(f + 4) (4 - *4) = ^(44 + 44 - «44 - ihh) 

= -{ij 5 hi 3 + ij%i 3 - 7 5 M 2 - 7 5 / / 2 ) 

= ~7 5 ^(44 + 44 - iloh ~ ihh) = -7 5 <?+e- • 

In a similar manner 



q-e+ = -7 5 9-e+ • 
9+e+ = l 5 q+e+ ■ 

From (jUJ) - dHJ) it follows that 

(1 + 7 5 )g + e_ = , 

(l + 7 5 )g_e + = 0, 
(l-7 5 )?+e + = 0, 
(1 -7 5 )g_e_ = . 



Hence 



(1 + 7 5 )e_e + «?_ = -2(1 + 7 5 )g_ , 
(l + 7 5 )e+e_g+ = -2(l + 7 5 )g+ , 
(l+7 5 )g-g+e- = 2(l + 7 5 )e_, 
(l+7 5 )g + g_e+ = 2(l + 7 5 )e+ . 



(1 - 7 5 )e+M_ = -2(1 - 7 5 )g_ , (34) 

(1 - 7 5 )e_e+g + = -2(1 - 7 5 )g + , (35) 

(l-7 5 )^_e_ = 2(l- 7 5 )e_, (36) 

(l-7 5 )^g + e + = 2(l- 7 5 )e + . (37) 
For example, for (pffij) we have 

(1 + 7 5 )e-e + g_ = 2(e_e+)(l + 7 5 )<?- - (1 + 7 5 )e+M- , 



where the second term on the right-hand side equals zero due to 

Furthermore, from the definition of the isotropic tetrad ([16)1 . (JT7j) it is clear that 

= 2<?_ , (38) 

4+9-<?+ = 2<?+ , (39) 

e_e+e_ = -2e_ , (40) 

e + e_e + = -2e + . (41) 

It is important, that any three operators of the form q\ or e a taken in succession in 
the presence of (1 ± 7 5 ) are either reduced to only one of them or equal zero. 

3 Recursion relations 

Now, using (J2IJ), <|2E|) - (SB) we have 

(1 + 7 5 )g_afo = (aq-)(l + J 5 )q-q+b - (ae+)(l + 7 5 )^e_fo 

= (og_)[(6gf + )(l + i 5 )q~q+q- - (6e+)(l + 7 5 )<?-<?+e~] 

- (ae+)[(6g_)(l + 7 5 )^e_g + - (6e_)(l + 7 5 )<?-e_e+] 

= (ag_)2[(6g + )(l + 7 5 )g_ - (6e + )(l + 7 5 )e_] 

- (ae + )2[-(fcg_)(l + 7 5 )e_ + (te_)(l + 7 5 )<7-] 

= 2[(ag_)(6g + ) - (ae + )(6e_)](l + 7 5 )<7- 

+ 2[(ae + )(6g_) - (ag_)(6e + )](l + 7 5 )e_ . (42) 
In the same way 

(1 + 7 5 )e_a6 = (og+)(l + 7 5 )<M-k - (ae_)(l + 7 5 )e_e+6 
= (og+)[(6g_)(l + 7 6 )M-£f - (be-) (I + 7 5 )M-e+] 
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- (ae_)[(6g+)(l + 7 5 )e-M- - (6e+)(l + 7 5 )e_e+e_] 

= (ag + )2[(6g_)(l + 7 5 )e_ - (6e_)(l + 7 5 )<7-] 
(ae_)2[-(6g + )(l + 7 5 )?- + (te+)(l + 7 5 )e-] 

= 2[(ag + )(6g_) - (ae_)(6e+)](l + 7 5 )e_ 

+ 2[(ae_)(6g + ) - (ag+)(6e_)](l + 7 5 )?- . (43) 
Denning the functions 

Fx (a, b) = 2[(aq^(bq + ) - (ae + )(6e_)] , (44) 

F 2 (a,6)=2[(ae + )(6 g _)-(a5_)(6e + )]=2 G ( ° M , (45) 



F 3 (a, b) = 2[(aq + ){bq_) - (ae_)(6e+)] , (46) 

F 4 (a,6)=2[(ae_)(6g + )-(a g+ )(te_)]=2 G( e a , (47) 

we finally get 

(1 + 7 5 )g_a6 = Fx(a, 6)(1 + 7 5 )g_ + F 2 (a, 6)(1 + 7 5 )e_ , (48) 

(1 + 7 5 )e_ab = F 3 (a, 6)(1 + 7 5 )e_ + F 4 (a, 6)(1 + 7 5 )g_ . (49) 
Making similar calculations, we obtain 

(1 + 7 5 )g + a6 = F 3 (a, 6)(1 + 7 5 )g+ + F 4 (a, 6)(1 + 7 5 )e+ , (50) 

(1 + 7 5 )e + a6 = F^a, 6)(1 + 7 5 )e + + F 2 (a, 6)(1 + 7 5 )g + . (51) 
In a similar manner, after defining the functions 

F 5 (a, 6) = 2[(aq-)(bq + ) - (ae_)(6e+)] , (52) 

F 6 (a,6)=2[(a e _)(6g_)-(a g _)(te_)]=2 G( e a ^) , (53) 

F 7 (a, 6) = 2[(ag + )(bg_) - (ae + )(6e_)] , (54) 

F 8 (a,6) = 2[(ae + )(6g+) - (« Z+ )(6e + )] = 2 G ( ° 6 ) (55) 

(1 - 7 5 )g_ab = F 5 (a, 6)(1 - 7 5 )g_ + F 6 (a, 6)(1 - 7 5 )e+ , (56) 

(1 - 7 5 )e + a6 = F 7 (a, 6)(1 - 7 5 )e + + F 8 (a, 6)(1 - 7 5 )g_ , (57) 
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we obtain 



(1 - r>)q + ab = F 7 (a, b)(l - 7 5 )g + + F 8 (a, 6)(1 - 7 5 )e_ , (58) 
(1 - 7 5 )e_a6 = F 5 (a, 6)(1 - 7 5 )e_ + F 6 (a, 6)(1 - 7 5 )g+ ■ (59) 
Note that the functions defined in (JH| - (J52j) - (|53j) are not independent: 

F 5 (o J 6) = J Fr(o,6) J 

F 6 (a,6) = F 2 *(a,&), 
F 7 (a,6) = F 3 *(a,6), 
F 8 (a,6) = F 4 *(a,6) . 
Furthermore, the following relations hold: 

F 3 (b,a)=F 1 (a,b) , 

F 1 (a,b)+F 3 (a,b)=2(ab) , 
Fi(a, 6)F 3 (a, 6) - F 2 (a, 6)F 4 (a, 6) = a 2 fe 2 

(see also Appendix B). 



4 A new method for calculation of traces 

We now use the formulae - (JoT?|) obtained in the previous section to calculate the 
traces of Dirac 7 -matrices as follows: 



Tr[(l + 7 5 )g_aia 2 a 3 • • • a 2n +i] 



= Fi(oi,a 2 )Tr 
+ F 2 (ai,a 2 )Tr 



Tr[(l + 7 5 )e_aia 2 a 3 • • • a 2n +i] = F 3 (ai,a 2 )Tr 

+ F 4 (ai,a 2 )Tr 

Tr[(l + 7 5 )g + aia 2 a 3 ---a 2n+ i] = F 3 (ai,a 2 )Tr 

+ F 4 (ai,a 2 )Tr 



Tr[(l + 7 5 )e+aia 2 a 3 • 



«2n+l] — Fi(tti, 

+ F 2 (ai,a 2 )Tr 



(1 + 7 5 )g_a 3 • ■■a 2n +i\ 
(1 + 7 5 )e_a 3 - • •a 2 „+i] 



a 2 ) Tr[(l + 7 5 )e+a 3 



(60) 



(1 + 7 5 )e_a 3 • • -a 2 n+i] 

(l + 7 5 )g_a 3 ---a 2n+1 ] , (61) 

(1 + 7 5 )g+a 3 • • -a 2n+ i] 

(l + 7 5 )e+a 3 ---a 2n+1 ] , (62) 



(1 + 7 5 )g+a 3 • • -a 2n+ i] , (63) 



and so forth. 
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As an example, let us calculate 

Tr[(l — 7 5 )aid 2 d 3 a 4 d 5 d 6 ] . 

We take decomposition (J2T|) written in the form 

(l + 7 5 )a 6 = (l + 7 5 )^(a 6 g + ) + (l + 7 5 )Q + (a 6 g_) 

- (1 + 7 5 )e-(a 6 e + ) - (1 + 7 5 )e+(a 6 e_) . (64) 

Using ([BU]) and (foTj) . we get 

Tr[(l + 7 5 )g_aid 2 d 3 d4a 5 ] 

= Fi(ai, a 2 ) Tr[(l + 7 5 )<2_d 3 a 4 d 5 ] + F 2 (ai, a 2 ) Tr[(l + 7 5 )e_d 3 a 4 d 5 ] 

= Fi(ai, a 2 ) {^1(03, a 4 ) Tr[(l + 7 5 )^a 5 ] + F 2 (a 3 , a 4 ) Tr[(l + 7 5 )e_a 5 ]} 

+ F 2 (ai, 02) {F 3 (a 3 , a 4 ) Tr[(l + 7 5 )e_a 5 ] + F 4 (a 3 , a 4 ) Tr[(l + 7 5 )9-fl5]} 

= [Fi(ai, a 2 )Fi(a 3 , a 4 ) + F 2 (ai, a 2 )F 4 (a 3 , a 4 )] • 4(a 5 g_) 
+ [Fi(ai, a 2 )F 2 (a 3 , a 4 ) + F 2 (a 1 , a 2 )F 3 (a 3 , a 4 )] • 4(a 5 e_) 

i.e. 

Tr[(l + 7 5 )g_aia 2 a 3 a 4 a 5 ](a 6 g + ) 

= 4[Fi(ai,a 2 )Fi(a 3 ,a 4 ) + F 2 (ai, a 2 )F 4 (a 3 , a 4 )] • (a 5 g_)(a 6 g+) 

+ 4[F 1 (ai,a 2 )F 2 (a 3 ,a 4 ) + F 2 (ai,a 2 )F 3 (a 3 ,a 4 )] • (a 5 e_)(a 6 g + ) . (65) 

Making similar calculations for the remaining three items in decomposition (|64|). we 
obtain 

Tr[(l + 7 5 )g + aia 2 a 3 o 4 a 5 ](a 6 g_) 
= 4[F 3 (ai,a 2 )F 3 (a 3 ,a 4 ) + F 4 (a 4 , a 2 )F 2 (a 3 , a 4 )] • (a 5 q + )(a 6 q_) 
+ 4[F 3 (ai,a 2 )F 4 (a 3 ,a 4 ) + F 4 (ai,a 2 )Fi(a 3 ,a 4 )] • (a 5 e+)(a 6 g_) , (66) 

- Tr[(l + 7 5 )e_aia 2 o 3 d 4 d 5 ](a 6 e + ) 

= -4[F 3 (ai,a 2 )F 3 (a 3 ,a 4 ) + F 4 (ai, a 2 )F 2 (a 3 , a 4 )] • (a 5 e_)(a 6 e+) 

- 4[F 3 (ai, a 2 )F 4 (a 3 , a 4 ) + F 4 (ai, a 2 )F 1 (a 3 , a 4 )] • (a 5 g_)(a 6 e+) , (67) 

- Tr[(l + 7 5 )e + did 2 a 3 d 4 d 5 ](a 6 e_) 

= -4[Fi(ai, a 2 )Fi(a 3 , a 4 ) + F 2 («i, a 2 )F 4 (a 3 , a 4 )] • (a 5 e+)(a 6 e_) 

- 4[Fi(ai, a 2 )F 2 (a 3 , a 4 ) + F 2 (ai, a 2 )F 3 (a 3 , a 4 )] • (a 5 g+)(a 6 e_) , (68) 
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and finally 

^ Tr[(l - rf)aia 2 a 3 a A a b a & ] 

= [F 1 (a 1 ,a 2 )F 1 (a 3 ,a 4 ) + F 2 (ai, a 2 )F 4 (a 3 , a 4 )]Fi(a 5 , a 6 ) 

+ [F^ai, a 2 )F 2 (a 3 , a 4 ) + F 2 (ai, a 2 )F 3 (a 3 , a 4 )]F 4 (a 5 , a 6 ) 

+ [F 3 (ai, a 2 )F 3 (a 3 , a 4 ) + F 4 (a 4 , a 2 )F 2 (a 3 , a 4 )]F 3 (a 5 , a 6 ) 

+ [F 3 (ai, a 2 )F 4 (a 3 , a 4 ) + F 4 (ai, a 2 )Fi(a 3 , a 4 )]F 2 (a 5 , a 6 ) , (69) 

that is we derive an expression consisting of eight terms only. 
An expression for 

Tr[(l + 7 5 )aia 2 a 3 a 4 d 5 a 6 ] 

may be obtained from 1)69)1 by replacing all the functions with the complex-conjugate 
ones. 

Having calculated 

Tr[(l - 7 5 )a 4 • • ■a 2n _ 1 a 2n ] , 
we may obtain an expression for 

Tr[(l - 7 5 )a 4 • • • a 2n _ia 2n a 2n+:L a 2n+2 ] 

from the previous one by the following replacement: 

Fi(a 2n -i, a 2n ) — > F\(a2n-i, «2n)Fi(a 2n+ i, a 2n+2 ) + F 2 (a 2n _i, a 2n )F 4 (a 2n+ i, a 2n+2 ) , (70) 
F 3 {a 2n -i-, o 2n ) — > F 3 (a 2n _i, a 2n )F 3 (a 2n+ i, a 2n+2 ) + F 4 (a 2n _i, a 2n )F 2 (a 2n+ i, a 2n+2 

), (71) 

F 2 (a 2n _i, a 2n ) — > Fi(a 2n _i, a 2n )F 2 (a 2n+ i, a 2n+2 ) + F 2 (a 2n _i, a 2n )F 3 (a 2n+ i, a 2n+2 

), (72) 

F 4 (a 2n _i, a 2n ) — > F 3 (a 2n _i, a 2n )F 4 (a 2n+ i, a 2n+2 ) + F 4 (a 2n _i, a 2n )Fi(a 2Tl+ i, a 2n+2 ) . (73) 
Actually, 

F 1 (a 2 n-i, a 2n ) = 2[(a 2n _ig_)(a 2n g+) - (a 2n _ie+)(a 2n e_)] 
= ^{ Tr[(l + 7 5 )g-a 2 „-i](a 2 „g+) - Tr[(l + 7 5 )e+a 2n „i](a 2 „e_)} 

^{ Tr t( 1 + 7 5 )?-a2n-ia2n«2n+i](o2n+2g+) - Tr[(l + 7 5 )e+a 2 „_ia 2 „a 2n+ i] (a 2n+2 e_) } 

= 2|[Fi(a 2 „_i,a 2n )(a 2n+ ig_) + F 2 (a 2n _i, a 2n )(a 2ri+ ie__)](a 2n+2 g + ) 

- [Fi(a 2 „_i,a 2 „)(a 2n+ ie + ) + F 2 (a 2n _i, a 2n )(a 2n+ ig+)](a 2n+2 e_)} 

— ^i(o 2n -i, a>2n)Fi(a 2n+ i, a 2n+2 ) + F 2 (a 2n _i, a 2n )F 4 (a 2n+ i, a 2n+2 ) , 

and the procedure for the remaining three functions is quite similar. 



Note that 



Tr[(l - 7 5 )aia 2 ] = 2(a 1 a 2 ) = a 2 ) + F 3 (a 1 ,a 2 ) 



(74) 



Thus, formulae (J74j) . (|7U|) - (|73|) provide a very simple method of trace calculation (see 
Appendix A). 

In numerical calculations, after fixing a reference system and choosing an orthonormal 
basis, for example, of the form 



# = (1, 0,0,0), %=(0, 1,0,0), # = (0, 0,1,0), Z£ = (0, 0,0,1) 



we get 



q± 



V2 



1, ±1,0,0) , 



4 



1 



(0,0,1, ±2), (oe±) 



and 



V2" 

Fi(a, 6) = (aft) - (a 6 2 



1 

i 



(o =F a*) 
a„ ± ia. 



a z & y ) , 



i(a y b z - a z b y ) , 



a^o) + i(a y b z 

F 3 (a, b) = (ab) + (a b x - a x b ) 

(a b y - a y b ) + i(a x b z - a z b x ) + (a x b y - a y b x ) + i(a b z - a z b ) 

(a b y - a y b ) + i(a x b z - a z b x ) - (a x b y - a y b x ) - i(a b z - a z b ) 

that is we derive a number of the expressions well-suited for the implementation of a 
simple yet efficient computer algorithm. 

Comparison of different methods for trace calculation is given in Tab. ^ 



F 2 (a,b) 
F A {a, b) 



(75) 
(76) 
(77) 
(78) 



Table 1: Number of terms in the expression for traces 



Number of 7-matrices, 2n 


2 


4 


6 


8 


10 


12 


Number of terms 
for algorithm (Jl]) 


Tr(aia 2 ■ ■ -a 2n ) 


1 


3 


15 


105 


945 


10395 


Tr(7 5 aia 2 • • • a 2 „) 





1 


15 


210 


3150 


51975 


Number of terms 
for algorithm (J2J) 


Tr 4 (aia 2 • • • a 2n ) 


1 


3 


15 


105 


693 


4383 


Tr 4 (7 5 aia 2 • • • a 2n ) 





1 


6 


33 


180 


1029 


Number of terms for 
new method (17011-11751) 


Tr[(l - 7 5 )aia 2 • • • a 2n ] 


2 


4 


8 


16 


32 


64 



In case that traces of the form 

Tr(7 P aia 2 • ■ ■ a 2n +i) Ti{Ybib 2 ■ ■ ■ b 2m+1 ) , 
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need to be calculated, one may use the consequence of formula (J2U|) 



1p x 7 



■p — 



e_ x e + — e + x e 



or sum over index p using the classical methods, e.g.: 



Tr(7 P aia 2 • ■ ■ a 2n +i) Tr(Yhb 2 ■ ■ ■ b 2m+1 ) 

= 2 Tr[(aia 2 • ■ • a 2n+1 + a 2n +i • • • a 2 ai)bib 2 ■ ■ ■ b 2m+1 ] . 



Finally, we would like to note that for traces of four, six and eight Dirac 7-matrices 
we have directly compared the expressions for the traces obtained both by the new and 
classical methods: 

- expressions for the new method were expanded by means of formulae (J75j) - (|75jl for 
functions F; 

- in expressions for traces obtained by the classical method the scalar products and 
contractions of the Levi-Civita tensors with 4- vectors were written in terms of vector 
components. 

The results obtained in different ways were identical. 

Note that some other examples of applications of isotropic tetrads to the calculations 
in high energy physics may be found in [7j. 



-Tr[(l - 7 5 )aia2a 3 a 4 ] = Fi(oi, a 2 )F 1 (a 3 , a 4 ) + F 2 (a 1 , a 2 )F^(a 3 , a 4 ) 

+ F 3 (oi, a 2 )F 3 (a^ a 4 ) + F 4 (ai, a 2 )F 2 (a 3 , a 4 ) ; (80) 

\ Tr[(l - 7 5 )aia2d 3 04a 5 a 6 ] 

= F x {a u 02)^(03, a 4 )Fi(a 5 , a 6 ) + F 2 (a 3 , a 4 )F 4 (a 5 , a 6 )] 

+ F 2 (ai,a 2 )[F 3 (a 3 ,a 4 )F 4: (a 5 ,a( i ) + F 4 (a 3 , a 4 )Fi(a 5 , a 6 )] 

+ F 3 (ai,a 2 )[F 3 (a 3 ,a 4 )F 3 (a 5 ,a 6 ) + F 4 (a 3 , a 4 )F 2 (a 5 , a 6 )] 

+ F 4 (a 1 ,a 2 )[F 1 (a 3 ,a 4 )F 2 (a 5 ,a 6 ) + F 2 (a 3 , a 4 )F 3 (a 5 , a 6 )] 

= [Fi(ai, a 2 )Fi(a 3 , a 4 ) + F 2 (a u a 2 )F 4 (a 3 , a 4 )]Fi(a 5 , a 6 ) 

+ [F 3 (ai, a 2 )F 4 (a 3 , a 4 ) + F 4 (oi, a 2 )F 1 (a 3 , a 4 )]F 2 (a 5 , a 6 ) 

+ [F 3 (ai, a 2 )F 3 (a 3 , a 4 ) + F 4 (a x , a 2 )F 2 (a 3 , a 4 )]F 3 (a 5 , a 6 ) 

+ [Fi(ai,a 2 )F 2 (a 3 ,a 4 ) + F 2 (ai,a 2 )F 3 (a 3 ,a 4 )]F 4 (a 5 ,a 6 ) ; (81) 



A Expressions for traces 



- Tr[(l - 7 5 )aia 2 ] = F^ai, a 2 ) + F 3 (ai, a 2 ) ; 



(79) 



1 
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forth. 



1 
2 



+ 



+ 



+ 



Tr[(l — 7 5 )oia 2 a^d 4 a5060 7 a 8 ] 

F 1 (a 1 ,a 2 )F 1 (a 3 ,a 4 ) + F 2 (ai, a 2 )F 4 (a 3 , a 4 
F ± (a 5 , a 6 )Fi(a 7 , a 8 ) + F 2 (a 5 , a 6 )F 4 (a 7 , a 8 

F 3 (oi, a 2 )F 4 (a 3 , a 4 ) + F 4 (ai, a 2 )Fi(a 3 , a 4 
Fi(a 5 , a 6 )F 2 (a 7 , a 8 ) + F 2 (a 5 , a 6 )F 3 (a 7 , a 8 
F 3 (ai, a 2 )F 3 (a 3 , a 4 ) + F 4 (ai, a 2 )F 2 (a 3 , a 4 
F 3 (a 5 , a 6 )F 3 (a 7 , a 8 ) + F 4 (a 5 , a 6 )F 2 (a 7 , a 8 

Fi(ai,a 2 )F 2 (a 3 ,a 4 ) + F 2 (ai, a 2 )F 3 (a 3 , a 4 
F 3 (a 5 , a 6 )F 4 (a 7 , a 8 ) + F 4 (a 5 , a 6 )F 1 (a 7 , a 8 



Tr[(l + 7 5 )o 1 o 2 • • • a 2n ] = (Tr[(l - 7 5 )o 1 o 2 • • • o 2n ])* 



Functions Fk(a, b) 



Fi(o,6) = 2[(og_)(6g+) - (ae+)(6e_)] = (06) + G 



F 4 (a, 6) 



o 6 



- G 



a 6 
^0 h 



i G 



a b 



e+ q. 
a b 



G 



h h 



= -lY[(l- 7 5 )g_e + a6] , 



2[(ae_)(6g + )-(a 0+ )(6e_)]=2 G 
iTr[(l- 7 5 )g + e_afe] . 



o 6 
e_ o 4 
o b 
h k 



+ 1 G 



= 1 Tr[(l - 7 5 )g + g-a&] = Tr[(l - 7 5 )^e + a6] , 

F 3 (a,b) = 2[(og + )(6g_) - (ae_)(6e + )] = (aft) - G -<G 

= 1 Tr[(l - 7 5 )g-g+afc] = Tr[(l - 7 5 )e + e„a6] , 

F 2 (o,6) = 2[(oe + )(6g_)-(ag_)(6e + )]=2G( e a + *) 



G 



+ i G 
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F 5 (a,b) = F*(a,b), F 6 (a, b) = F*(a, 6), F 7 (a, b) = F 3 *(a, b), F 8 (a, b) = F 4 *(a, b). (88) 
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